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Compensated Convexity, Multiscale Medial Axis Maps 
and Sharp Regularity of the Squared Distance Function 

Kewei Zhang* Elaine Crooks^^and Antonio Orlando^ 


Abstract 

In this paper we introduce a new stable mathematical model for locating and measuring 
the medial axis of geometric objects, called the quadratic multiscale medial axis map of scale A, 
and provide a sharp regularity result for the squared-distance function to any closed non-empty 
subset K of K". Our results exploit properties of the function C'pdist^(-; K)) obtained by ap¬ 
plying the quadratic lower compensated convex transform of parameter A |59j to dist^(-; K), the 
Euclidean squared-distance function to K. Using a quantitative estimate for the tight approx¬ 
imation of dist^(-; K) by (7pdist^(-; K)), we prove the C^’^-regularity of dist^(-; K) outside a 
neighbourhood of the closure of the medial axis Mk of K, which can be viewed as a weak Lusin- 
type theorem for dist^(-; K), and give an asymptotic expansion formula for C'pdist^(-; K)) in 
terms of the scaled squared distance transform to the set and to the convex hull of the set 
of points that realize the minimum distance to K. The multiscale medial axis map, denoted 
by M\{-; K), is a family of non-negative functions, parametrized by A > 0, whose limit as 
A —>■ oo exists and is called the multiscale medial axis landscape map, Moo(-; K). We show that 
Moo(-; K) is strictly positive on the medial axis Mk and zero elsewhere. We give conditions 
that ensure K) keeps a constant height along the parts of Mk generated by two-point 

subsets with the value of the height dependent on the scale of the distance between the generat¬ 
ing points, thus providing a hierarchy of heights (hence, the word ’multiscale’) between different 
parts of Mk that enables subsets of Mk to be selected by simple thresholding. Asymptoti¬ 
cally, further understanding of the multiscale effect is provided by our exact representation of 
Moo{-', K). Moreover, given a compact subset K of M", while it is well known that Mk is not 
Hausdorff stable, we prove that in contrast, M\{-] K) is stable under the Hausdorff distance, 
and deduce implications for the localization of the stable parts of Mk- Explicitly calculated 
prototype examples of medial axis maps are also presented and used to illustrate the theoretical 
findings. 
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1 Introduction 

The medial axis of an object is a geometric structure that was introduced by Blum m as a means 
of providing a compact representation of a shape. Initially defined as the set of the shock points of 
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a grass fire lit on the boundary and allowed to propagate uniformly inside the object, closely related 
definitions of skeleton m and cut-locus [53] have since been proposed, and have served for the 
study of its topological properties laEaiiiiiiiiEii, its stability pdlEI] and for the development 
of fast and efficient algorithms for its computation m im dn Eg ssj. Applications of the medial 
axis are ample in scope and nature, ranging from computer vision to image analysis, from mesh 
generation to computer aided design. We refer to [50] and the references therein for applications 
and accounts of some recent theoretical developments. 

An inherent drawback of the medial axis is, however, its sensitivity to boundary details, in the 
sense that small perturbations of the object (with respect to the Hausdorff distance) can produce 
huge variations of the corresponding medial axis. This observation has prompted a large body of 
research that has roughly followed two lines, both aimed at the definition of some stable modification 
of the medial axis: one consists of reducing the complexity of the medial axis by pruning the less 
important parts of the domain |3g, the other considers the definition of filter conditions that 
identify subsets of the medial axis which are stable to perturbations of the sets and retain some of 
its topological properties, for instance, homotopy equivalence with the object. Within this second 
line of research, we mention, among others, the A—medial axis introduced in (20], the 0—homotopy 
preserving medial axis introduced initially in |3T| and subsequently modified in [52] to ensure the 
homotopy equivalence, and the power crust method [S|. The A-medial axis and the 0-homotopy 
preserving medial axis are explicitly defined as subsets of the medial axis, being collections of those 
points of the medial axis that meet some geometrical criteria. Such criteria are expressed in terms 
of a bound either on the distance to the boundary of the object or on the separation angle 9 (see 
Definition 3. 16] below), respectively. The power crust, in contrast, and also the algorithm discussed 
in pg |26| , provide a continuous approximation of the medial axis constructed using a subset of the 
vertices, called poles, of the Voronoi diagram of a finite point sample of the object boundary. In all 
such works, the stable modifications are sought by identifying directly points of the medial axis or 
of an approximation of it. The excellent survey paper m contains a thorough discussion of such 
approaches and of the related stability issues. 

We adopt in this paper a fundamentally distinct strategy which, if compared with the works 
mentioned above, represents an indirect approach relying on the use of the compensated convex 
transforms |59|. The theory of compensated convex transforms has been introduced and applied 
in the calculus of variations for finding the quasiconvex envelope of a function [Sg Eg EZl EH] 
and for finding tight smooth approximations of the maximum function and the squared-distance 
function [60] . Compensated convex transforms, however, also provide a natural and stable global 
method to extract geometric singularities, such as ridges, valleys and edges, from a given function 
by manipulating its ‘landscape’ [621163] . and it is in this way that the transforms, in particular the 
lower compensated convex transform (hereafter, called also the lower transform), will be used in this 
paper. Whether one applies the lower compensated convex transform or the upper compensated 
convex transform depends on the type of geometric singularities to be extracted. The works [621 
63] present a systematic study on the use of these basic transforms to extract singularities from 
the graph of functions in general, or from the characteristic functions of compact sets, whereas 
the patent application m contains various applications including our method for extracting the 
multiscale medial axis map. The key properties that are exploited to highlight and/or to design 
a specific singularity are: the tight approximation of the compensated transforms, their regularity 
and the manner in which they respond to the type of curvature. More specifically, [62] focuses on 
the basic use of these transforms to detect ridges, valleys and saddle points of graph of functions, 
whereas m presents the design of a transform which is capable of filtering out the ‘regular points’ 
and the ‘regular directions’ on manifolds. 

The application of the lower transform to study the medial axis of a set is motivated funda- 
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mentally by the identification of the medial axis with the singularity set of the distance function 
m Lemma 8.5.12] and by the geometric structure of this set [3l HSl 03]. On the other hand, the 
distance function, its regularity and its geometric structure, are well-studied both in geometric 
measure theory m and in the theory of partial differential equations [H [291 132 EZj. If the set 
K is a smooth compact submanifold of M”, there are many local regularity results of the distance 
function near K gni Ea EH ED, whereas, for a general bounded open set 0, some results by Al- 
bano la imply that the distance function dist(-; 11'^) is locally in \ in the sense that if 
X £ VL \ Mqc, there is a 5 > 0 such that dist(-; G J)). 

In the following, however, it is more convenient to refer to the squared-distance function and 
use the identification of the singular set of the distance function with the set of points where the 
squared distance function fails to be locally Here, we just note that the advantage of referring 
to the squared distance function rather than to the distance function has also been realized in 
other contexts, such as, in the study of the motion of surfaces by its mean curvature represented 
by manifolds with codimension greater than one [22 El- We refer to |2 for a detailed study on 
the properties of the squared distance function and on its applications in the geometric evolution 
problems. 

Using properties of the lower transform, we apply the lower compensated convex transform to 
the Euclidean squared-distance function which gives a smooth tight approximation outside a 

neighbourhood of the closure of the medial axis (see Theorem 3.3), and define our multiscale medial 
axis map as a scaled difference between the squared-distance function and its lower transform. From 
the property of the tight approximation of the lower transform of the squared-distance function, 
we also deduce a sharp C^’^-regularity result (see Corollary 3.8 and Example 3.10) of the squared- 
distance function outside a neighbourhood of the closure of the medial axis of K, which can be 
viewed as a weak Lusin type theorem for the squared-distance function and extend regularity results 
of the squared-distance function to any closed non-empty subset of M”. This result also offers an 
instance of application of the compensated convex transform to obtain a fine result of geometric 
measure theory and is, somehow, related to the behaviour of semiconcave functions (see m and 
c) below). We observe that, in general, the regularity of dist^(-; K) cannot be better 
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Remark 

than even for a compact convex set K where Mk = 0. A simple example is given by the 
square K = [0,1] x [0,1] C M^. In this case, it can be easily verified that dist^(x; K) is globally 
but not C2. 

The application of the lower compensated convex transform of scale A to the squared-distance 
function produces a continuous function in that remains strictly positive on the medial axis and 
tends to zero outside of it as a positive parameter A becomes very large (see Proposition 3.20). We 
will, in fact, characterize the limit of the multiscale medial axis map of scale A as A approaches to 
inhnity (see Theorem 3.23) and refer to this geometric structure as the quadratic multiscale medial 
axis landscape map of K. The values of this map are well separated, in the sense that they are 
zero outside the medial axis and remain strictly positive on it. Furthermore, we will give conditions 
(see Proposition 3.18 and Section]^ that ensure that the multiscale medial axis map of scale A 
actually keeps a constant height along the parts of the medial axis generated by two-point subsets, 
with the value of the height dependent on the distance between the two generating points. Such 
values can, therefore, be used to define a hierarchy between different parts of the medial axis and 
we can thus select the relevant parts through simple thresholding, that is, by taking suplevel sets 
of the multiscale medial axis map. To reflect this property, we use the word ’’multiscale”. For each 
branch of the medial axis, the multiscale medial axis map automatically dehnes a scale associated 
with it. In other words, a given branch has a strength which depends on some geometric features 
of the part of the set that generates that branch. 

Given a closed non-empty subset K of M”, we will also prove that, despite the medial axis of K 
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not being Hausdorff stable, the quadratic multiscale medial axis map, is indeed Hausdorff stable 
(see Theorem 4.3). It follows that the graph of the medial axis map carries more information than 
the medial axis itself, which allows the definition of a hierarchy between the parts of the medial 
axis and the selection of the relevant ones through simple thresholding, that is, by taking suplevel 
sets of the medial axis map. In this manner, it is possible to choose the main parts that reflect 
genuine geometric features of the object and remove minor ones generated by noise. 

In conclusion, we observe that while our method seems to share similarities with those based 
on the extraction of ridges of the distance transform [9l dU |39l 051 [5l], (that require, however, 
an a-priori definition of ridge, based usually on an approximation of the derivative of the distance 
transform), the method we propose is, in fact, substantially different from such approaches, given 
that we obtain a neighbourhood of the singularities as the difference between the squared-distance 
transform and its smooth tight approximation. In this manner, as mentioned above, we provide an 
indirect definition of the singularity, which does not require any derivative approximation or any 
differentiability assumption. 

After this brief introduction, the next section will introduce the relevant notation and recall 
basic results in convex analysis and lower compensated convex transforms. Section contains the 
definition of the multiscale medial axis map, and some of its principal properties, such as the tight 
approximation of the lower compensated transform to the squared-distance transform (see Theorem 


3.3) and as an application, we deduce a sharp regularity result of the squared-distance function to 
any non-empty closed subset of M” (see Corollary |3.8[ ). Section [^presents the Hausdorff stability of 
the multiscale medial axis map whereas Section discusses some mathematical prototype models 
of explicitly calculated medial axis maps for a simple four point set to some more complicated three 
dimensional objects. Finally, Section concludes the paper with the proofs of the main results. 


2 Notation, Basic Definitions and Preliminary Results 

Throughout the paper M" denotes the n-dimensional Euclidean space, and |x| and x-y the standard 
Euclidean norm and inner product, respectively, for x, y € In some cases, we will also make use 
of the notation (x, y) to denote the point of given by xei + ^262 -f .. • -f j/nCn, where {ei,..., en} 
is an orthornormal basis of (x, y 2 , ■ ■ ■, yn) £ IK"' and y = 7/262 . -f y-aen- Given a non-empty 

subset K of M"', denotes the complement of K in M’'^, i.e. \K, K its closure and co[K] 

the convex hull of K, that is, the smallest (with respect to inclusion) convex set that contains the 
set K. For x G M” and r > 0, i?(x; r) indicates the open ball with center x and radius r whereas 
S{x] r) denotes the sphere with center x and radius r and is the boundary of B{x', r). The distance 
transform of a non-empty set K C M" is the function that, at any point x G associates the 
distance of x to K, which is defined as inf{|x — y\, y G K} and is denoted as dist(x; K). We use 
the notation Df to denote the derivative of /. 

Across the current literature, there is no uniform definition of the medial axis, with its meaning 
changing from one author to another. What the medial axis is for one, becomes the skeleton for 
another, and in some cases subtle differences are present, especially in the continuum case, where 
the closure of such sets is considered. In this paper we adopt the definition given by Lieutier in 
im, but it is reformulated here to include a non-empty closed set K C M"' with K 7 ^ M"' as well as 
a non-empty bounded open set Q. 

Definition 2.1. For a given non-empty closed set K C with K 7 ^ MF, we define the medial 
axis Mk of K as the set of points x G M" \ AT such that x G Mk if and only if there are at least two 
different points 7/1, 7/2 £ K , satisfying dist(x; K) = \x — yi\ = \x — 7 / 2 |- For a non-empty bounded 
open set Q. C M", the medial axis of Q is defined by Mq := n Mqq. 
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Remark 2.2. (a) The definition of medial axis for a bounded open set 0 is equivalent to that of 

the closed set since the definition of Mk implies that Mqc = Mqq, and hence Mq = Mqc. 

(b) Other frequently used notions are those of the skeleton of K, denoted as skl{K), and the cut 
locus of a manifold, denoted as cl{K), which applies to the more general case of Riemannian 
geometry. Here, for a non-empty closed subset K of the Euclidean space M”, we define the 
skeleton of K to be the set of the centers of the maximal (with respect to inclusion) open balls 
contained in K^, whereas the cut locus of K is taken to be the set of the cut locus of the 
points of the boundary of K in M”", where the cut locus of p in dK is the set of points in the 
manifold where the geodesics starting at p stop being minimizing. It can then be shown that 
cl{K) = Mk. As a result, the notions of medial axis, skeleton and cut locus are related but 
are not the same; for example, in general, [2^ 

Mk C skl{K) cMfi. (2.1) 


(c) Our definition of the medial axis Mk is, in particular, consistent with our main convergence 
result (Theorem 3.g5| ) which recovers Mk as the set on which the limit of the medial axis map 
M\(-; K) (Definition 3.1) as X tends to infinity is strictly positive. 


Next we collect definitions and results from convex analysis for functions / taking finite values, 
i.e. for / : M” i—)■ M, which will be used in this paper, and refer to [361 Wf\ for details and proofs. 


Given a function / ; i—?• M bounded below, the convex envelope co[/] is the largest convex 

function not greater than /. We will often make use of the following characterization. 


Proposition 2.3. Let f : M"’ i—)> M 6e coercive in the sense that f(x)/\x\ —)> oo as |x| —>■ oo, and 
xq G M". Then 

(i) The value co [/] (xq) of the convex envelope of f at xq G M"" is given by 


(ii) 


n+l 


n+1 


n+l 


CO [/] (xo) = _ inf < ^ Xif{xi) : X] ^ X] ^ *0’ A* > 0, x* G M” > . (2.2) 

*-i,...,n+i y j 

If, in addition, f is lower semicontinuous, the infimum is reached by some (A*, x*) for i = 
1, 2,..., n + 1 with (x*, f{x*)) ’s lying in the intersection of a supporting plane of the epigraph 
of f, epi(/), and epi(/). 

The value co [/] (xq), for f taking only finite values, can also be obtained as follows: 

CO [/] (xo) = sup{t'(xo) : i affine and l{y) < f{y) for all y G M”} (2-3) 

with the sup attained by an affine function i* G Aff(M’^). 


We then introduce the following definition |13) . which is needed in Lemma 2.5 


Definition 2.4. Assume xq G M”. We say that a function f : M" i-+ M is upper-semidifferentiable 
at Xo if there exists a G M" such that 


lim sup 
y^O 


f{x + y) - f{x) -a-y 

■ ly\ 


< 0 , 
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The following lemma, concerning the existence and properties of an optimal affine function, will 
be needed for the proofs of Proposition 2.14 and Theorem |3.3[ 


Lemma 2.5. Suppose f : M"" i—?• M is continuous, upper-semidifferentiable, coercive in the sense 
that lim|3,|_,.oo /(x)/|x| = + 00 , and co[/] G //co[/](0) < /(O), then there is an affine 

function i{x) = 2a ■ x + b and distinct points xi,X 2 , ■ ■ ■ ,Xk G M"' and Ai > 0, A 2 > 0, ..., Afc > 0 
satisfying 2 < k < n + I, Yli=i Yli=i — 0 Xi xj if 1 < i j < k, such 

that 

(i) i{x) < f{x) for all x G ; 

(ii) i{xi) = f{xi) for i = 1,2,... ,k; 

(in) 2a = Di{xi) = Df{xi) for i = 1,2,... ,k; (2.4) 

[iv) b = £{0) = co[f]{0) ■, 

(u) 2a = T>£(0) = T»co[/](0). 


The quadratic lower compensated convex transform, introduced in [59], will play a pivotal role 
in the definition of our multiscale medial axis map. We next recall its definition and some of its 
properties, and refer to [59l ESI E3| for details and proofs. 

Definition 2.6. Let f : M”’ t—>■ M 6e a lower semicontinuous function satisfying 

f{x)>-Ai\x\‘^-A2 (2.5) 

for some constants Ai, A 2 >11. The (quadratic) lower compensated convex transform (lower trans¬ 
form for short) for f with scale A > j4i is defined for x G M"' by 

Ci{f){x) = co[f + A| • p](x) - A|xp (2.6) 

If f is bounded below, we may set A > 0. 


Remark 2.7. (a) The requirement of the lower semicontinuity of f is to guarantee that C\{f){x) —)> 

f{x) as A —>■ 00 for all x G MT, since otherwise, the lower transform will converge to the lower 
semicontinuous envelope of f. 


(b) From (2.6) it also follows that C'^(/)(x) is the envelope of all the quadratic functions with 
fixed quadratic term A|xp that are less than or equal to f, that is, for x G M” fUR Eq. (l.f)] 


C'^(/)(x) = sup {—A|xp + t'(x) : —A|yp + £(y) </(y) for all y G M” and f affine}. (2.7) 


(c) Recalling from m that a function f : M is called c—semiconvex if, for some constant 

c > 0, the function /(x) + c/2|xp is convex, we observe that the lower compensated convex 
transform for f with scale X, C^if) is a 2X-semiconvex function. In fact, Cgf) represents the 
2X-semiconvex envelope of f. We sometimes use such a property to extend some properties 
of semiconvex functions to C\{f). 

(d) To gain further geometric insight into the lower compensated convex transform defined by 
( |2.6| ), in Fiyure we display the steps of the construction of C\^{f) for /(x) = dist^(x; K) 
with K = { — 1, 1} and A = 2. The graph of the augmented function / + A|xp is displayed 
in Figureg^b) along with f, whereas Figureg^c) shows the convex envelope of the augmented 
function. Figure S'*) displays finally the graph of C'}(/) which is compared with that of f. 
Note that the convex envelope of the augmented function is different from f + A|xp only in 
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a neighbourhood of the singular point 0 of f, so that, when then we subtract the weight, the 
final effect is a smoothing of f only in such neighbourhood. This simple example, along with 
the ones discussed in Section enables one also to understand the role of the parameter 
A and our meaning of scale. The parameter A acts as a scale parameter in the sense that 
it controls the curvature of the lower compensated convex transform in the neighbourhood of 
the singularity of the function and allows the extraction of the singularity with a value which 
gives somehow a measure of its strength. Also one may observe the so-called ’tightness’ of the 
lower compensated convex transform approximation of the original function from below (see 
Proposition 2.10), which agrees with the original function except in the neighbourhood near 
the singular point. 



Figure 1: Steps illustrating the construction of the lower compensated convex transform of f{x) = 
dist^(x; K) with K = —1, 1. (a) Graph of the function /; (6) Graph of the augmented function 
/ + A| • p with A = 2; (c) Graph of the convex envelope of / + A| • p compared to that of / + A| • p; 
{d) Graph of C\^{f) compared to that of /. 


The following properties of C\^{f) will also be used. 

Proposition 2.8. Given f : M"" i—)> M that satisfies (2.5), then for all Ai < \ < t < oo, we have 


Ci(/)(x) < C|(/)(x) /orxGM", 


( 2 . 8 ) 
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and, for A > ^41 


(2.9) 


C{{f){x) < f{x) for X G M” . 

Proposition 2.9. If f < g in M"' and satisfy ( |2.5[ ), then 

C[{f){x) < C{{g){x) for x G M” and A > maxj^ij, Ai^g} . (2.10) 


The transform C\^{f) realizes a ‘tight’ approximation of the function /, in the following sense 
(see [59l Theorem 2.3(in)]). 

Proposition 2.10. Let f G C^’^{B{xo] r)), with B{xo‘, r) the open ball of center xq and radius 
r > 0. Then for sufficiently large X > 0, we have that /(xq) = C'^(/)(xo). 


Such a property motivates the definition of the mnltiscale ridge transform which was introduced 
in [62] to extract ridges of general fnnctions and shown to be invariant with respect to translation. 
This multiscale ridge transform will be used in Section]^ to define the multiscale medial axis map 
(see Definition 3.1). 


Definition 2.11. Given A > 0, the ridge transform of scale X, for a given function f : 




satisfying (2.5), is defined as: 


Rxifix))-.= f{x)-C[{f){x) 


x G 


( 2 . 11 ) 


We now present some regularity properties of C'((_(dist^(-; K)), which will be exploited to analyze 
the behaviour of the multiscale medial axis map. We recall first the following result given in |59( 
Lemma 4.3]. 

Lemma 2.12. Suppose f : M” — >■ M fs convex and such that |/(x)| < cjxp + ci for x G M"" with 
c, Cl > 0. Assume A > 0 and define /a = A| • p — /. Then for X > c, 

co[fx]ix + y)-co[fx]{x)-Dco[fx]ix)-y< X\y\^ (2.12) 


for x,y €MT. 

The next proposition improves a result in |59[ Theorem 3.1]. 

Proposition 2.13. Suppose K C M” is a non-empty closed set. Then for A > 0, 

Ci(dist2(.; K)) G ^^’^(M"). 

Furthermore, the Lipschitz constant of the gradient iAC'(^(dist^(-, K)) is at most 2max{l, A}. 

The next property is a useful inequality for the derivative of the lower transform, IlC'^(dist^(-; K)). 
Proposition 2.14. Suppose K C M” is a non-empty closed set. Then 

\DCiidist\-, K)){x)\^ < 4Ci(dist2(.; K)){x) (2.13) 


and equality holds in (2.13) if and only i/C'^(dist^(-; K)){x) = dist^(x; K). 





Remark 2.15. Given that u = dist^(-; K) is known to be a viscosity solution of the Hamilton- 
Jacobi equation EEI 

\Du\^ = Au, (2.14) 


Proposition 2.13 and Proposition 2.14 together imply that the lower transform C^(dist^(-; K)) is a 
subsolution of the Hamilton-Jacobi equation (2.14). 


We next introduce the sets K{x) and K 2 ^\{x), which will be used to gain insight into the 
geometric structure of C'^(dist^(-; K)). 

Definition 2.16. Let K C M"' he a non-empty closed set. For any x G M"', let r(x} = dist(x; K). 
We then define the following sets: 


and for A > 0, 


iir(x) = |x + y G M” : x + y G dK and \y\ 


1^2,\{x) := ^ G M” : x + y G dK and 

’ L 1 + A 


= r(x)| 

\y\ = ^(®)} • 


(2.15) 

(2.16) 


Remark 2.17. If x ^ K, the set K{x) is the set of points of dK that realize the distance of x to 
K. Note also that it follows from (2.15) that K{x) = B{x; r{x)) n K, so in particular, K{x) is 
compact, and if x ^ K, B{x; r(x)} C K‘^. 


The following result, obtained in the proof of m Theorem 3.7], gives an explicit expression of 
the lower transform of dist^(-; K{x)), the squared distance to the set K{x), and will be used to 
produce a bound on the multiscale medial axis map (see Theorem 3.15[i)). 


Proposition 2.18. Let K C M"" be a non-empty closed set and Mk the medial axis of K. Assume 


X G Mk and denote by K{x) and K 2 ^\{x) the sets defined by (2.15) and by (2.16), respectively. 
Then, for all y gW^, 


A 


CA(dist (•; K{x))){y) = (1 + A)dist (y; co[il' 2 ,A(a;)]) + (x) - X\y - x 


(2.17) 


where co[K 2 ^x{x)] is the convex hull of K 2 ^x{x). 

We will also need, for the proof of Theorem |3.23[ the following explicitly calculated formula of 
the lower transform for compact sets contained in a sphere 5(0; r) = {x G jxj = r} centred at 
0 G M"' with radius r > 0. The formula is easy to derive following similar calculations to those in 
the proof of |60l Theorem 1], or of |59l Theorem 5.1]. 

Proposition 2.19. Let K C 5(0; r) be a non-empty compact set. Then for every x G M”, 


Ci(dist2(.; K))ix) = 


Ar^ 

1 + A 


+ (1 + A)dist^ 


x: CO 


K 

1 + A 


— AJx] 


(2.18) 


where co[ii'/(l + A)] is the convex hull of K/{1 + A) = {x/(l + A), x G K}. 


We will invoke the following technical lemma several times (see Lemma 3.2 in |59j ). 
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Lemma 2.20. Assume p > 0. Let K = p) = {x G \x\ > p} he the eomplement of the 

open ball -B(0; p) with center the point 0 and radius p, then 


C'i(dist2(.; K)){x) 


A 


-p^ — A|x| 


1 + A' 
dist^(x; K), 


x\ < 

|x| > 


1 + A’ 
P 

1 +A ■ 


(2.19) 


In the next lemma, which generalizes slightly [591 Lemma 3.3], we give the expression of the 
lower transform of the squared distance to a set of two points. The two points, without loss of 
generality, are assumed to lie along a basis vector of specifically, along the basis vector ei G 
This lemma will be used extensively when we investigate the behaviour of the multiscale medial 
axis map with respect to perturbations of the boundary of K. 

Lemma 2.21. Assume n >2 and let {ei, ..., Cn} he an orthonormal basis of the Euclidean space 
M"". Let K = {—aei, aei}, where a > 0. We write y = 622/2 + • • • + e„ 2 /n £ and represent the 

point xei + y £ as the pair (x, y), which therefore denotes the point (x, 2/2, • ■ ■, yn) G M”. Then 
for every A > 0, we have 


Ci(dist^(-; K)){x,y) 


- Ax^ + |2/|^ 

dist^((x, 2 /); K) 


< 

> 


a 

iTa’ 

a 

1 + A ■ 


In particular, 


C'i(dist^(-; K)){0,y) = + 1^1^ < + |2/|^ = dist‘^{{0,y); K). 


( 2 . 20 ) 


( 2 . 21 ) 


We conclude this section with the definition of 5—neighbourhood of a set, of Hausdorff distance 
between two sets |6], and of e—sample of a set [TO] . 

Definition 2.22. Given a non-empty subset E ofW^ and 6 > 0, we define the 5-neighbourhood E^ 
of E by 

E^ = {x gM."' : dist(x; E) < <5}. 

Note that E^ is an open subset of M"'. 

Definition 2.23. Let E, F be non-empty subsets ofMA. The Hausdorff distance between E and F 
is defined in m by 

dist^(.B, F) = inf |5 > 0 : T C and E C (2.22) 

This definition is also equivalent to saying that 


distu{E,F) 


max 


< supdist(x; F), supdist(x; E) > . 

1 x^E TCP 1 


xGF 


Definition 2.24. Let K he a compact subset of ME. A sample S of the boundary of K is a finite 
set of points of the boundary of K, i.e. S C dK and #(5) G N where ff{S) denotes the cardinality 
of S. An e—sample of dK is a sample whose Hausdorff distance to dK is less than e, that is, 
dist'p(S', dK) < e. 

A uniform e—sample S of dK is an e—sample of dK such that 


dist^(S, SAT) < ediam(Ar), 


(2.23) 
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where the diameter of K, diam(iir), is defined as 

diam(iir) = sup \x — y\ . 

x,y&K 

3 The Multiscale Medial Axis Map 

In this section, we define the quadratic multiscale medial axis map Mx(-;K), characterize some of 
its properties, and establish its relation to the medial axis Mk- As a by-product, we also infer sharp 
regularity results for the squared distance function dist^(-;-A), which are of independent interest. 

Definition 3.1. Let K C be a non-empty closed set. The qnadratic mnltiscale medial axis 

map of K (medial axis map for short) with scale X > 0 is defined for x G M"' by 

Ma(x; K) := (1 + A)iiA(dist2(-; K)){x) = (1 + A) (dist2(x; A)) - Ci(dist2(-; A)))(x)) . (3.1) 

For a bounded open set Ll C with boundary dLl, we define the qnadratic multiscale medial 
axis map of Ll with scale A > 0 as 

Ma(x; n) := Ma(x; dQ) x £ Q. (3.2) 


Remark 3.2. (a) The convergence of the lower transform to the original function as A —?• oo 

yields that limA^oo .RA(dist^(-; A))(x) = 0, implying that the values of the ridge transform 
can be very small when X > 0 is large. To make the height of our medial axis map on the 
medial axis bounded away from zero, we thus need to scale the ridge transform. The factor 
(1 -|- A) turns out to be the “right” scaling factor, as will be justified in Theorem 3.15 below, 
where it will be shown that on the medial axis Mk, the medial axis map Mx{x] A) is bounded 
both above and below by quantities independent of X. 


(6) The quadratic multiscale medial axis map can also be seen as a morphological operator ra 
equal to the scaled top-hat transform of the squared distance transform with quadratic structur¬ 
ing function. Letting f{x) = dist^(x; A) and bx{x) = —X\x\‘^, it can he shown that the lower 
transform corresponds to the grayscale opening operator with quadratic structuring function 

m; i-e-, 

c{{f) = {febx)®bx, 


and thus. 


Ma(-; A) = (l + A) 


f - (febx)® bx 


Notwithstanding such an interpretation, it is convenient to view Definition \3.1\ in terms of 
the lower compensated convex transform. The exploitation of properties of such transforms 
permits a relatively easy evaluation of the geometrical properties of Mx{-', A) and also permits 
an easy numerical realization of Mx{-] A). This relies on the availability of numerical schemes 
for computing the lower transform of a given function, which entails the availability of schemes 
to compute the convex envelope of a function. We refer to for the algorithmic and 
implementation details of the schemes for realizing the lower transform of a function. 


We begin with a key quantitative estimate of the tight approximation of the squared distance 
function dist^(-; A) by its lower transform C'^(dist^(-; A)). This result not only underpins our 


11 





study of the role of Mx in characterizing the medial axis Mk, but also yields improved locality 
and regularity properties of C'^(dist^(-; K)) and dist^(-; K) respectively (see Corollaries 


3.6 


3.1 


and 


3.13) which are of interest in their own right. 


Theorem 3.3. Let K C M" he a non-empty closed set and denote by Mx the medial axis of K. 
Suppose Mk 7 ^ 0 , a > 0, assume x G M" \ Mx, and let Mx{x; K) he the multiscale medial axis 
map of K with scale A. If 


^ ^ dist(x; K) 

~ dist(x; Mx) ’ 


(3.3) 


then 

dist^(x; K) = C'^(dist^(-; K)){x), 


(3.4) 


and consequently, 


Mx{x] K) = 0. 


(3.5) 


Remark 3.4. Note that Mx = 0 if and only if K is convex (see, for example, Theorem 2.21]), 
in which case dist^(-; K) is convex, and therefore equals C')^(dist^(-; K)) in for all A > 0. 


Now assume A > 0 and introduce the set 


Vx,K = {x G M"' : Adist(x; Mx) < dist(x; K)} . (3-6) 

Clearly Mx C Vx^X: so this defines a “neighbourhood” of the medial axis Mx of K (note that it 
is possible that Mx n ill / 0 , so VA,ic is not necessarily a neighbourhood in the strict sense), and 
Vx,x is a closed set. Moreover, as A > 0 increases, Vx,x describes a family of shrinking sets such 
that 

[^Vx,x=^, (3.7) 

A>0 

and if we take the support of the multiscale medial axis map. Theorem |3.3| yields that 


spprt(MA(-; K)) cVx,x ■ 


(3.8) 


With the help of (3.8), we can show the following result that characterizes Mx iu terms of 
spprt(MA(-; K)). 

Corollary 3.5. Suppose K C M” is a non-empty closed set and Mx 7 ^ 0. Then 


spprt(MA(-; 7^)) = Mk. (3.9) 

0<A<+oo 


An important consequence of Theorem 3.3 is the following locality property of the lower trans¬ 
form of the squared distance function. This result is also of independent interest, in particular 
because it quantihes the size of neighbourhood needed to evaluate C')^(dist^(-; 7C)), and also be¬ 
cause it will be exploited in the proofs of Theorems 3.15 and 3.15 to establish results characterizing 
the properties of Mx{-;K). 


Corollary 3.6. (Locality Property) Suppose K C M” is a non-empty closed set. Then for every 
xq G M”, 

Ci(dist2(.; K)){xo) = C 05 (,,o;^(,,o))[dist 2 (-; K) + A|(-) - xo|^](xo), (3.10) 
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where 


(3.11) 


r(xo) := -dist(xo; K). 


Remark 3.7. (a) CoTollaTy \3.^ improves the result in f38^ . where the radius of the ball for the 

locality property is r{xo) = 6dist(xo; K)/^/\ for A > 2. 

(b) In IMf, it was also established that xq G M"" is a stationary point 0 /C'(^(dist^(-; K)) if and 
only if xq G co[K(xo)]. We will see that the ‘only if’ part of this result is a consequence of 


arguments from the proof of Theorem 3.3: see Remark 6.1 


Theorem 3.3 can also be combined with Proposition 2.13 to yield a regularity property of the 


distance transform, which can be viewed as a weak version of the Lusin theorem for the squared- 
distance function. 

Corollary 3.8. Assume A > 0. Let K CW"' be a non-empty closed set and V\^k the neighbourhood 


of Mk defined by (3.6). Then 


dist2(.; K) GC^’l(M'^\PA,i^). 
Furthermore, for all x, y £ \ V\^k 

|Zldist^(y, K) — iAdist^(z, K)\ 


\y - ^1 


< 2 max{l, A} . 


(3.12) 


Remark 3.9. It follows from Corollary 3.8 and Stepanov’s Theorem 
dist^(-; K) exists almost everywhere in M"" \ Mk- 


that the Hessian of 


Both estimate (3.3) in Theorem 3.3 and estimate (3.12) for the Lipschitz constant in Corollary 


3.8 (when A > 1) are, in fact, sharp, as the following example shows. 


Example 3.10. Consider K = (—1, l)'^ = M \ (—1, 1) C M. Then Mk = {0} and for A > 0, 

1 


Ci(dist2(.; K)){x) = 


A , 2 , , 

— Ax , |x| < 


A + 1 
dist^(x; K), 
0 , 


A + 1’ 


1 


A + 1 
|x| > 1. 


< kl < 1, 


As a result, dist^(x; K) = C'^(dist^(-; K)){x) if and only if |x| > 1/(A + 1). Moreover, for this 
example we have that dist(x; Mk) = \x\, whereas dist(x; K) = 1 — |x| if \x\ < 1 and dist(x; K) = 0 

is given by 


if \x\ > 1. Therefore (3.3) holds if \x\ > 1. If 0 < |x| < 1, 

A|x| >1 — |x|<t4>|x| >1/(1 +A), 


which shows that estimate (3.3) is sharp for K = (—1,1)'^. 


Furthermore, for x± := ±1/(A + 1), 

Zldist^(x±, K) = +4A/(1 + A), |x+ — x_| = 2/(1 + A) 
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so that 


|Z)dist^(x+, K) — L>dist^(x_, K)\ 


= 2A. 


|rE+ — x_| 

Thus we can also conclude that estimate (3.12) of Corollary \3.^ is sharp when A > 1. 

We cover next the case of a bounded open set C M”, giving a precise statement about equality 


of the medial axis maps Ma(x; and M\[x] dQ), followed by a modification of Corollary 3.8 
Proposition 3.11. Suppose H C MT is a non-empty bounded open set and let A > 0. Then 

C'^(dist^(-; dn)){x) = C'^(dist^(-; Tl^)){x), x G (3.13) 

and consequently, 


M\{x] 511) = M\{x; QS), x G H. 


(3.14) 


Remark 3.12. (a) Property (3.14) of the medial axis map is important in many practical situ¬ 

ations. For example, in image processing, the objects H of which we wish to find the medial 
axis might be defined by taking a threshold from a greyscale image, that is, as a suplevel set 
of the image function. The object is then represented by a binary image rather than by its 
boundary. Therefore, in this case, it might be more convenient for us to compute numerically 
the medial axis map Mx{x] H'^) rather than M\{x\ 511). 

{h) It is worth noting the different qualitative behaviour of the convex envelope and the compen¬ 


sated convex transform that appears in (|3.13|). 
always positive in 


For A > 0, the left hand side of (3.13) is 


By setting A = 0, 


\ n whereas the right hand side equals zero in 
the left hand side of (3.13) reduces to the the convex envelope o/dist^(x; 511), which vanishes 


in the convex hull co[ll] of the closure of II, whereas the right hand side of (3.13) gives the 
convex envelope o/dist^(x; H'^), which is identically zero in M"'. 


For a bounded open set H, Corollary 3.8 modifies as follows. 
Corollary 3.13. Let n C M"" be a bounded non-empty open set. Then 

dist2(.; H^) GC^’i(ll\WA,n), 

where 

Wx,n := jic G H, dist(x;MQc) < ^ 

and diam(ll) is the diameter of H. Furthermore, for all x, y £ Il \ Wx,n, 

\DdistHy, !!“) - Odist^(^ ^ 


\y - ^1 


(3.15) 

(3.16) 

(3.17) 


A consequence of Corollary 3.13 is that outside any neighbourhood of dist^(-; H'^) is a (7^4 


function. However, we also notice that Mk can have positive n-dimensional Lebesgue measure. 


Therefore the measure of Wx,n might not be small even when A > 0 is large. Corollary 3.8 and 
Corollary 3.13 also demonstrate that the lower transform can be viewed as a (7^4 extension of the 
squared-distance function from the set Vxk^ on which dist^(-; AT) = (7^(dist^(-; iT)), to M"' and from 
fl \ Wx,n to 11, respectively. 


14 


















Theorem 3.3 showed that if x ^ Mk, then M\{x\K) = 0 when A is sufficiently large. We now 


further explore the relationship between the medial axis map and the medial axis Mj^, 

both establishing A-independent positive upper and lower bounds on M\{x] K) whenever x G Mk, 
and fully characterizing the limit of M\{-‘,K) as A —)• oo. The following geometric structure will 
play a key role in both Theorem 3.15|and Theorem 3.23 


Definition 3.14. Let K C 


be a non-empty closed set and for x € denote by K{x) the set 
defined by (2.15), that is, K{x) = B{x] r{x))r\K, and denote by co[/C(a:)] the convex hull of K{x). 
The quadratic multiscale medial axis landscape map of K is defined for x G ML by 


Moo(x; K) := dist^(x; K) — dist^(x; co[K{x)]) 


(3.18) 


It is straightforward to see that Moo(x; K) = 0 if x 0 Mk but Moo(x; K) > 0 for all x G Mk- 
Indeed, if x 0 Mk, then there exists y & K such that K{x) = {y} = co[it'(x)] and dist^(x; K) = |x — 
yp, thus Moo(x; K) = 0. On the other hand, if x G Mk, then there exist distinct yi, y 2 G K{x), so 
since Ki+M g co[K{x)], we have Moo(x; K) > ^jx—yip + ^|x—y 2 p —|x — > 0. 

The next result establishes key bounds on Mx{-,K). 


Theorem 3.15. Let K C ML be a non-empty closed set, and denote by Mk the 
and by Moo(x; K) the quadratic multiscale medial axis landscape map defined by 


medial axis of K 


(i) For every A > 0 and every x G Mk, 


Moo(x; K) < Ma(x; K) < dist2(x; K). 


(3.19) 


(ii) For every A > 0 and for every x G 

0 < M\{x; K) < dist^(x; K) . 


(3.20) 


The lower bound in (3.19) can be expressed in terms of the separation angle which has been 
used, for instance, in [52], for a local geometrical characterization of the medial axis. 


Definition 3.16. Let K C be a non-empty closed set and denote by Mk the medial axis of K. 
For X G Mk, let yi, y 2 G K{x) and denote by Z[yi — x,y 2 — x] the angle between the two non-zero 
vectors yi — x and y 2 — x, taken between 0 and vr, i.e. Z[yi — x,y 2 — x] = cos“^ 
then define the separation angle Ox for x G Mk as follows: 


Ox 


max 



- x,y 2 - x]. 


yi, 02 G 



(3.21) 


Remark 3.17. Recall from Remark 2.11 that K{x) is compact and hence the supremum of the set 
{Z[yi — X, y 2 — x], yi, y 2 G K{x)} is realized by a pair of distinct points of K{x). 


Proposition 3.18. Let K C M”' be a non-empty closed set, and denote by Mk the medial axis of 
K. Then for every A > 0 and x G Mk, 

sin^(0a:/2)dist^(x; K) < Ma(x; K) < dist^(x; K). (3.22) 
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Remark 3.19. (a) For a limit point x G Mk \ Mk, we have either x £ K or Ox = ^, and the 

estimate ( |332l ) is valid also in this case. 


(b) Since the bounds in (3.22) are independent of X, we have that for x G Mk, 


s\ti{9x/2) dist^(x; K) < liminf Ma(x; K) < limsup Mx{x; K) < dist^(x; K). (3.23) 

A^+oo A^+oo 


(c) Since for x G Mk, there are at least two different points yi, y 2 £ dB{x; r{x)) H K, it follows 
that Ox > 0, and hence, comparing with (3.22), we have Mx{x; K) > 0 for x £ Mk- 


The next result gives the limit behaviour of Mx{x', K) as A —)• cx) for x 0 Mk- 

Proposition 3.20. Let K C M"' be a non-empty closed set and denote by Mk the medial axis 
of K. Assume A > 0 and denote by Mx{x] K) the medial axis map of K of scale A. Then for 
X G M” \ Mk, 

lim Mx{x] K) = 0 . (3.24) 

A—^~1“00 


Remark 3.21. Theorem 3.3 clearly implies (3.24) special case that x 0 Mk. But the set 


Mk \ {Mk U K) may not be empty and w e do not know whether Mx{x, K) > 0 for all X > 0 for 
X £ Mk \ {Mk U K). So Proposition 3.20 is needed for the proof of the general characterization of 
lim;,^oo Mx{-', K) that will be given in Theorem 3.23. 


Remark 3.22. By Proposition 3.18 and Proposition 3.2C\ we have 


lim inf Mx{x; K) 
A—^“)~oo 


> 0, if X £ Mk , 
= 0, if X ^ Mk . 


This implies that the limit of Mx{x; K) can extract exactly the medial axis of K. If we apply a 
slightly weaker scaling to the ridge transform, say (1 + x)“ for 0 < a < 1, and define M“(x; K) = 
(1 + x)“iiA(dist^(-; K)){x), we have 


{ +CX) , if x £ Mk , 

= 0 , if X ^ Mk , 

that is, M“(x; K) approaches the indicator function of Mk ]36f as X becomes large. 

We can now characterize the limit of Ma(x; K) as A —)> +oo for all x G M"". 
Theorem 3.23. Suppose K C M”' be a non-empty closed set. Then for every x £ 

lim Mx{x; K) = Moo{x] K). 

A—^ “1“CXD 


(3.25) 


From Theorem |3.3t Corollary 3.8 and Theorem 3.23 it follows that when A > 0 is increasing, 
the support of Mx{-', K) is contained in a shrinking neighbourhood of Mk and approaches the 
multiscale medial axis landscape map Moo(x; K). The numerical advantage of studying Mx{-; K) 
as an approximation of the multiscale medial axis landscape map Moo{-; K) is that it relies only 
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on the computation of the lower compensated convex transform of the squared distance transform, 
whose construction is local by virtue of Corollary 3.6, whereas the computation of Moo{-; K) is 
difficult because we need to evaluate the convex hull co[ifr(x)]). 


Remark 3.24. A further consequence of Theorem \3.23[ is that for every fixed x G M" and for every 
non-empty closed set K C M"’, the family of lower transforms A i—)• C^(dist^(-; K)){x) is ‘differen¬ 
tiable’ at infinity. If we let e = 1/X, g{e] x) = C'{y^(dist^(-; K)){x), and g{0-, x) = lime_^o x) = 
dist^(x; K), then 


ff(e; x) -g(0; x) 
e-s-0+ e 


-Moo(x; K), 


so we have the asymptotic expansion 


C^(dist^(-; K)){x) = dist^(x; K) 


Moo(x; K) 

l + A 



(3.26) 


when A —)• oo. 


Remark 3.25. In general, x i—)• dist^(x; co[iir(x)]) is not continuous in M” as x approaches the 
medial axis Mk- But from Theorem \3.23[ we can show that for every x G M", 

^lim |(1 + A)C'j^(dist^(-; K)){x) — Adist^(x; ifr)| = dist^(x; co[it'(x)]) (3.27) 

using the equality 

dist^(x; K) — Mx{x-, K) = {1 -\- A)C'^(dist^(-; K)){x) — Adist^(x; K). (3.28) 

For large A > 0, ( |337l ) can be viewed as an approximation o/dist^(x; co[it'(x)]) by continuous 
functions. The function dist(x; co[77(x)]) has been used, for instance, for surface reconstruction 
when K is finite IWf . While it is difficult in general to calculate dist^(x; co[ifr(x)]) directly, 

we see from (3.27) that the numerical computation o/C')^(dist^(-; K)){x), whose evaluation involves 
only local convex envelope calculations because of Corollary \3.(^ offers an easy approximation of 
dist^(x; co[ifr(x)]). 


We conclude this section by observing briefly that, based on the estimates of Theorem 3.15 and 
Proposition 3.18[ it is reasonable to define an alternative medial axis map by taking the square 
root of Ma(x; K). 

Definition 3.26. We define the mnltiscale medial axis map of linear growth (linear medial 
axis map for short) by 

M{(x; K) := a/Ma(x; K) (3.29) 

for X GMF. 


From Proposition 3.18, we obtain that the height of this linear medial axis map is ‘proportional’ 
to the distance function itself; that is, for A > 0, we have 


and 


sm{9x/2) dist(x; K) < M\{x] K) < dist(x; K) for x G Mk 
0 < M\{x; K) < dist(x; K) for x G M”. 


(3.30) 

(3.31) 
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Note that the linear medial axis map M|(x; K) is different from a definition based on the 
lower compensated convex transform for the distance function dist(x; K) itself, i.e. based on 
i?A(dist(-; K)). Of course, we can define such maps using the p-distance function dist^(rc; K) for 
any 1 < p < oo. But in this paper, we focus mainly on the medial axis map defined 

using the squared-distance function, i.e., for p = 2, in which case the geometry of M\{x] K) is easy 
to control. For instance, as we will see in the next section, M\(x; K) has the same height along 
the parts of the medial axis generated by two points. This is a key property when one looks for 
approximate medial axes by applying the Voronoi diagram method of finite e-samples. 


4 Hausdorff Stability 


Quantifying the instability of the medial axis is of fundamental importance for both theory and 
computation. This aspect becomes more and more relevant in practice nowadays, given that point 
clouds are increasingly being used for geometric modeling over a wide range of applications. More¬ 
over, there are computational approaches, such as the Voronoi diagram method, which search for 
a continuous approximation of the medial axis of a shape starting from subsets of the Voronoi 
diagram of a sample of the shape boundary. The presence of noise on the boundary, and/or the 
discrete character of samples of the boundary shape thus call for methods that permit the control 
of the parts of the medial axis which are not stable. In this section, we will discuss how this aspect 
is tackled by the multiscale medial axis map. In the first part of the section, we examine the values 
of M\{-] K) when the distance of the point to the boundary of the set is achieved by two points, 
whereas in the second part we discuss the Hausdorff stability of M\(-] K). 


Proposition 4.1. Assume n>2 and let {ei, ...,en} be an orthonormal basis of the Euclidean 
space Let K = {—aei, aei}, where a > 0. We write y = 62^2 + • • • + enVn £ and 

represent the point xei + y ^ as the pair {x,y) G M x which therefore denotes the point 

{x,y 2 , ■ ■ ■, yn) £ Then for every X > 0, we have 


Mx{{x,y); K) 


(1 -F A)^ 

0 , 



x\ < 
x\ > 


a 

iTa’ 

a 

1 + A’ 


(4.1) 


Remark 4.2. The medial axis map M\{{x,y); K) reaches its maximum on the medial axis of K, at 
the point (0, y) = ^262 + - • • + ynen £ attaining the value Ma((0, y); K) = o?. Note that a > 0 

is half the distance between the two points —aei and aei of K. Another important observation is 
that Mx{{x,y); K) is a function of the x-variable only, and does not change its height along the 
y direction. Therefore, on branches of the medial axis generated by two points, the height remains 
the same. If a is small (equivalently, the two points in K are close to each other), the values of 
Mx{{x,y); K) will be uniformly small. 

We next give the Hausdorff stability property of the multiscale medial axis map, followed by 
some comments on implications of this property for the localization of the medial axis of a domain. 


Theorem 4.3. Assume A > 0. Let K, L d'ST be non-empty compact sets. Then as L ^ K under 
the Hausdorff distance, Mxf; L) —)■ Mx{-; K) uniformly in every fixed bounded set in M”. More 
precisely, if we let fi ;= dist'i^(ii', L) he the Hausdorff distance between K and L, then for x G M”, 


Ci(dist^(-; K)){x) - L)){x) 


< h 


X', K) /i)^ -|- 1 -|- /i^ 


(4.2) 


18 



and 


M\{x] K) — M\{x; L) < ^{1 + A) ^(dist(x; K) + + 2dist(x; K) + 2^ + 1^ . (4.3) 


Remark 4.4. (a) From Theorem 4-3 we also conclude that for any compact sets Ki, K 2 

Ma(x; Ki) - M\{x; K 2 ) < min |/i(l + A) (^(dist(x; Ki) + p.)'^ + 2dist(x; Ki) + 2/i + | 


which shows that the medial axis map is uniformly continuous on compact sets with respect 
to the Hausdorff metric. 


{b) While the medial axis of K is not a stable structure with respect to the Hausdorff distance, its 
medial axis map M\{x\ K) is by contrast a stable structure. This result complies with (4.3) 
which shows that as A becomes large, the bound in (4.3) becomes large. 


As an immediate consequence of Theorem 4.3 we have the following result, which relates the 
medial axis map of the boundary of a domain Q with that of its e-samples K^. 


Corollary 4.5. Assume A > 0. Let Ll C M” be a bounded open set with diameter di'\air){Ll). Suppose 
Kf: C M” is a compact set such that dist'^(9n, AT^) < e. Then for x G M"' 


C'A(dist^(-; dQ)){x) — C'A(dist^(-; Kf)){x) < e^(diam(n) -|- e)^ -|- 1 -|- e^ , 


and 


Mx{x; dQ) — M\{x; Kfj < e(l -|- A)^(diam(r2) -|- e)^ -|- 2diam(n) -|- 2e -|- 1^ , 


for all X € Q. 


(4.4) 

(4.5) 


Remark 4.6. (a) If we consider an e-sample AT^ of dH, that is, a discrete set of points such that 

dist'^(5n, ATe) < e, Corollary \4. ^ yields a simple criteria that permits the suppression of those 
parts of the Voronoi diagram of AT^ that are not related in the limit, as e —)■ 0, to the stable 
parts of the medial axis of LI. 


(b) Since the medial axis of AT^ is the Voronoi diagram of Kf., if denotes the set of all the 
vertices of the Voronoi diagram Vor{Kf) of K^,, and is the subset of V^ formed by the 
poles ofVor{K,:) introduced in (i.e. those vertices ofVor{Kf) that converge to the medial 
axis of Ll as the sample density approaches infinity), then as a result of Proposition 3.20. for 
X > 0, we conclude that 


lim M\{xe, KA = 0 for G 14 \ 

e->0+ 


(4.6) 


Since as e ^ 0-|-, ATg —)• dLl, and knowing that P^ —)• 0 \wi, then on the vertices of 

Vor(A’e) that do not tend to Mq, M\{Xf:', Kf) must approach zero in the limit because of 
Proposition \3.26 . As a result, in the context of the methods of approximating the medial axis 
starting from the Voronoi diagram of a sample (such as those described in l^\2^\2^\5Df), the 
use of the multiscale medial axis map offers an alternative and much easier tool to construct 
continuous approximations to the medial axis with guaranteed convergence as e ^ 0-|-. 


With the aim of giving insight into the implications of the Hausdorff stability of M\{x', dLl) 
and Corollary |4.5t we display in Figure the graph of the multiscale medial axis map of a non- 
convex domain Ll and of an e-sample A’e of its boundary. Inspection of the graph of Ma(x; dLl) and 
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M\{x] Ke), displayed in Figure j^o) and Figure [^6), reveals that both functions take comparable 
values along the main branches of Mq. Also, M\[x] takes small values along the secondary 
branches, generated by the sampling of the boundary of O. These values can therefore be hltered 
out by simple thresholding so that a stable approximation of the medial axis of can be computed. 
This can be appreciated by looking at Figure [^d), which displays a suplevel set of M\{x\ K^) that 
appears to be a reasonable approximation of the support of Ma(x; dO) shown in Figure j^c). 



Figure 2: Multiscale Medial Axis Map of a nonconvex domain 0 and of an e-sample of its 
boundary, (a) Graph of Ma(-; 0) for A = 5; (6) Graph of M\{-] K^); (c) Support of M\{-] Vt)] (cl) 
Suplevel set of Mx{x] K^) for a threshold equal to 0.15 max{M;,(rE; Ale)}. 

a;eK2 


5 Examples of Exact Medial Axis Maps and Their Supports 

In this section we illustrate the behaviours of our multiscale medial axis map for some 2d geometric 
objects K, for which it is possible to obtain an explicit analytical expression for M\{x; K). Thanks 
to the translation and the partial rotation invariance property of the convex envelope [621 Proposi¬ 
tion 2.3, 2.10], it is then possible to derive an explicit analytical expression for M\[x] K) in the case 
that AT is a 3d-solid obtained by, for instance, rotations or translations of the models considered in 
this section. For the sake of conciseness, we leave the derivations to interested readers. 
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Though the derivation here is limited only to 2d geometric models, these models retain, never¬ 
theless, their basic geometric features, because they are able to show that M\{x; K) can, in fact, 
provide an accurate and stable way to find Mk, the medial axis of K, and represents likewise an 
effective tool to analyze the geometry and structure of Mk- We will also see how it is possible 
to select either the main stable parts of or to locate its fine parts by using suplevel sets of 
Ma(x; K). 

Example 5.1. We consider the case of a four-point set K C defined as follows. Let 6, e > 0 
with e G (0, 1), set c = eb. Define then K = {(6, c), {b,—c), {—b,c), (— 6, —c)}. For this set, we 
have 

dist2((x, y); K) = (|x| - hf + (|y| - cf , (5.1) 

and, after some lengthy calculations based on the construction of affine functions, we can show that 
the lower transform can be expressed as follows 

Ci(dist^(-; K)){x,y) = {l + X)g{x,y) + + c^) - Xix"^ + y"^) (5.2) 


where the auxiliary function g = g{x,y) is a continuous piecewise quadratic function defined as 
follows 


g{x,y) = < 


0 


X — 


\y\ - 


X — 


1 + A 


1 + A 


1 + A 


+ |y|- 


1 + A 


if |x| < 
if |x| > 
if |x| < 
if |x| > 


1 + A’ 
b 

TTa’ 

b 

1 + A, 
b 

TTa’ 


\y\< 
\y\< 
|y|> 
\y\ > 


1 + A’ 
c 

iTa’ 

c 

iTa’ 

c 

1 + A 


(5.3) 


The multiscale medial axis map M\{x\ K) is then computed using the definition (3.1). In particular, 
for this example, after some algebraic rearrangements, it is possible to show that since all four points 
in K lie on a circle centerd at the origin, the medial axis map of K can be expressed as 


G 


(5.4) 


M\{{x,y); K) = {1 + A)^(^dist^((x, y); K/{1 + A)) - dist^((x,y); co{K)/{l + A))^ (x,y) 

where for Z C M” and a G M \ {0}, we use the notation Zja to denote the set {w G M” : w = 
zja for z G Z}. By a closer inspection of (5.4), we can make then the following observations: 

(i) The support of M\{-; K) is 


spprt(MA(-; K)) 


|(x,y) G : |x| < 


2b 

1 + A 


or |y| < 


2eb 1 

iTa/ ■ 


The ‘thickness’ of the support for the main branch y-axis is, therefore, 26/(1 + A) while that 
for the minor branch x-axis is 2e6/(l + A). 

{ii) The height of the medial axis map along the main branch y-axis when \y\ > e6/(l + A) is b^ 
while the value along x-axis when |x| > 6/(1 + A) is e‘^b‘^. 

{Hi) At the only Voronoi vertex Q of the Voronoi diagram of K, the value is Mx{0; K) = 6^(l + e^). 
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Figure^ displays the graph of Mx{x; K) as given by (5.4) for different values of X and for the set 
K defined by b = 2 and e = 0.5. For each value of X, we can easily verify the presence of two scales 
in M\{x; K): a strong one which is reflected by the values of M\{x; K) along the y—axis generated 
by the two-point set {{—b, c), {b, c)} and a weak one captured by the value of M\{x] K) along the 
x—axis generated by the two-point set {{b, —c), (—6, —c)}. In agreement with our theoretical results, 
we also verify that the support of the continuous function M\{x‘, K) contains the medial axis of K, 
given by the Voronoi diagram of K in this case, and such support shrinks to Mk as X increases. 



Height 



Hei^ 


{d) 



Width 


Height 




Figure 3: Graph of Ma(-; K) for the four-point set K = {{b,c), (6, —c), {—b,c), (— 6, —c)} with 
6 = 2 and c = e6 = 1, for different values of A: (a) A = 0.25; (6) A = 4; (c) A = 8; (d) A = 15. 


Example 5.2. In this example, we consider first the case of the open set Hg = (—r, oo) x (—r, r) 
with r > 0, whose results will be used to construct the multiscale medial axis map of a rectangular 
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domain. By inspection, we can easily infer that 


dist^((x,y); dQg) = < 


i\y\-rf 

{\y\-r? 

{\x\-rf 

{\x\-rf 


if X > O', 

if — r < X < 0, |y| > |x|; 
if -r <x <0, \y\ < |x|; 
if X < -r, \y\ < r; 


(5.5) 


_ (|x| - r)^ + (|y| - r)^ if x <-r, |y|>r, 

whereas the lower transform, obtained after lengthy ealculations based on the construction of affine 
functions, is given, for {x,y) G Qs by 


Ci(dist^(-; dns)){x,y) = 
g{x + r,y + r) 
g{x + r,r-y) 


if -r<x<0, -r<y, x + y < 

r 

if — r < X < 0, —r< —y, x — y<—- 


X 


1 + A 


1 + A 


- r 


— A(x^ + if X <0, — 


1 + A 


<x + y, - 


1 + A 


1 + A’ 
<x-y; 


-Ay' 


if x>0, \y\ < 
if X > 0, 


1 + A’ 

< |y|; 


(5.6) 


1 + A 

{\y\ - rf if X> -r, \y\ > r; 

(|x| — r)' if X < —r, \y\ < r; 

^ (|x| - r)'+ (|y| - r)' if x <-r, \y\ > r , 

where the auxiliary function g = g{x,y) is a continuous piecewise quadratic function defined as 

follows 


g{x,y) = 

^ 2 

X 


if X < 


A 


a; > 0, y> 0; 


1 + A’ 

X A 

if y < y>o, x>0; 


1 + A 




*/ <y < y > 0 , a; > 0 ; 


1 + A 
if X <0, y > 0 


A 


y */ X > 0, y < 0; 

^ x' + y' if X <0, y < 0 , 


The multiscale medial axis map of Vlg, Mx{{x,y)', Ttg), is obtained by applying definition (3.1) 
and by taking into account (5.5) and (|5.6|). By exploiting properties of the lower transform with 
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respect to symmetry and translation of axis, we can then easily obtain the analytical expression for 
the multiscale medial axis map of a rectangular domain. If, for instance, we consider the open 
bounded set Q = 2 ) ’ ^ 2 ) ^ difficult to show that 


Mx{{x,y)] Q) = 


^x{{x +-,y)-,Ils) ifx<Q, 
Mx{{-x+^-,y)-,VLs) i/x>0. 


(5.8) 


Figure^^a) displays the support of Mx{{x,y); Q) which is a neighbourhood of the medial axis, 
whereas Figure^^b) depicts the graph of M\{{x,y); bl). For the points {x,y) G Mq with Ox = tt, 
it follows that M\{{x,y)] Q) = dist^((x,y); dQ), so implying that in this sense, the upper bound in 
(3.22) is sharp. 



Figure 4: Multiscale medial axis map of O = {— {r + ^) ,r + x (—r,r) for A = 10: (o) Graph 
of Mx{{x,y); fl) for A = 10; (6) Domain Q displayed along with the support of Mx{{x,y); D). 


Example 5.3. We consider now the oval shaped domain D C made by the union of two semi¬ 
circles with center at the points (—r/2, 0) and [rj^, 0) and radius r jl, respectively, and the rectangle 
{—rl2,rl2) x (—r/2, r/2). For this domain, it is not difficult to verify that 


dist^((x,y); dbl) = < 


(|y| - rf if |x| < 



(5.9) 


whereas the lower transform, obtained after some lengthy calculations based on the construction of 
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affine functions, is given, for {x, y) G by 


C'i(dist^(-; dn)){x,y) = 



i\y\-rf 



r 

X + -<0; 
r 

X + -<0; 
r 

i-2>o; 

X - - > 0; 
2 “ 


(5.10) 


Figure^a) displays the graph of Mx{{x,y); Q) obtained by applying definition (3.1) where we 
account for (5.9) and (5.10), whereas Figure^b) shows the support of M\[{x,y)] Q) along with the 
domain 11. 



Figure 5: Multiscale medial axis map for the oval shaped domain: (a) Graph of the medial axis 
map, M\{{x,y); fl), for A = 10 and of the squared distance function, dist^((x,y); 511) whose zero 
level set gives the boundary of the domain. Only the restriction to y > 0 is displayed, given that 
M\{{x,y)] n) < dist^((a:, y); 5H); (6) The oval shaped domain H displayed with the support of its 
medial axis map M\{{x,y)\ H). 


In the following example we describe the behaviour of M\(x; K) for the case of a discrete set 
K, sampled from a connected set, and evaluate the structure of M\{x] K) as the sample density 
approaches infinity. 

Example 5.4. We consider the geometric model of the uniform sampling of two parallel lines at 
distance b to each other. The points are taken equally spaced over each line at distance c = eb with 
e G (0,1) measuring the sampling density, in the sense that as e —>■ 0, the sampling density on the 
two lines tends to infinity. The sampling of the two parallel lines is defined so that the discrete 
points are aligned along the x— and y—axis as displayed in Figure^a). For such a sample K, we 
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can then use the results obtained in Example 5.1 for the four-point set which we refer to as K^p. It 
is not difficult to show that, for {x,y) G and i G Z such that \y + ie6| < eb, 


Mx{{x,y)-, K) = Mx{{x,y + ieb)] K^p ), 


(5.11) 


where Mx{{x,y + ieb); K^p) is the multiscale medial axis map for the four-point set discussed in the 
Example \5.1\ 



Figure 6: (a) Geometric model of a set K representing a uniform e—sample of two parallel lines. 

Graph of the multiscale medial axis map for different sample density: (6) 6 = 2; e = 0.5; A = 10; 
(c) 6 = 2; e = 0.1; A = 10. 


Figure and^c) display the graph of Mx{{x,y); K) for 6 = 2, A = 10 and for different 
different sampling density e, e = 0.5 and e = 0.1, respectively. 

The comparison of the graphs of Mx{{x, y)] K) for the two different sample densities of the same 
object, shows how the value of Mx{{x,y); K) along the minor branehes of the medial axis of the 
discrete set K attenuates os e —)> 0. Consistently with the finding obtained for the four-point set, we 
have that the value of Mx{{x,y); K) along the minor branches is proportional to e^. It then follows 
that by setting a threshold not lower than such a value, we can single out the stable part of Mk, 
which provides an approximation for (and in this case, is in fact coincident with) the medial axis 
of the two parallel lines. 

In this last example, we analyze a model of perturbations of the boundary domain represented 
by staircase-like piecewise affine curves. This effect is very common, for instance, in digital images 
and is the source of unrealistic medial axis branches, which are usually not desirable. Gommon 
practice in this case is to perform a boundary smoothing prior to any image processing operation. 
We will show that this is not needed with the multiscale medial axis map. The fine structure of the 
medial axis corresponding to the irregularities of the boundary is indeed captured by the multiscale 
medial axis map and can be filtered out. We will verify this statement on a prototype model of 
this boundary domain perturbation, by showing that the height of the medial axis map on such 
branches can be very small if the stair like effect is small. 

Example 5.5. Assume c > 0 and let us consider the set K = {{x,y) £ : x>0, y>0, x + y< 

c} n {{x,y) G M^, X -I y > c} displayed in Figure^ which is used as a prototype of one single step 
perturbation of a boundary domain. 
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Figure 7: Geometric model of one stair-like boundary perturbation, as in a digitized domain. 


The squared distance function dist^((x, y); K^)) to the complement of K and the lower transform 
C'j^(dist^(-; K^)){x,y) are then given, respectively, by 



if X < y, 0 < X < c, 0 < y < c. 


if X > y, 0 < X < c, 0 < y < c. 

+ {y- cf 

ify>c, 0 <y -X <c, 

(x - cf + 

if X > c, 0 < X — y < c, 

lix + y-cf 

if X + y > c, y — X > c, 


or X + y > c, X — y > c, 

0 

otherwise ; 


(5.12) 
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and 


C'i(dist2(-; K^)){x,y) = 
2/A 


0<a:<c, 0 < y < c, 


y 

A(l + A) 


if y < , 0<x<c, 0 <y < c, 


1 + A’ 


x + y 


— A X — 


1 + 2A 

Ax (1 + A)x 

if <y < 


{l + X){x + y) 

1 + 2A 


-Xly- 


(l + A)(x + y) 


1 + 2A 


1 + A 


0<x<c, 0<y<c, x + y< 


(1 + 2A)c 


1 


A(x + y-c)2; 
x^ + (y - cf 
(x - cf + 


if X + y > c, y — X > c, or x + y>c, x — y > c, 


ify>c, 


< y — X < c, 


1 + A 

Q 

if X > c, -- < X — y < c, 

j _ , i^x- ^ ’ 


Ac2 
1 + A 


- X{y - x)^ + (x + y - c)' 
«/a: + y > (1 + 2A) 


1 + A’ 1 + A 


> y — X > — 


1 + A’ 


2 I 2 

X +y 


«/a;<0, y < 0, 
otherwise. 


(5.13) 


The medial axis map M\{{x, y); K^) is then obtained from using ( |5.12[ ) and ( |5.13[). The graph 
of M\{{x, y); K'^) for A = 9 and step size c = 1 is shown in Figure^^a), whereas Figure^b) displays 
its support. By inspeeting the graph of Mx{{x,y); K^), we observe that after an initial inerease near 
the corner tip, M\{{x,y); K‘^) keeps a constant value along Mk, with this value proportional to the 
square of the step size. It is not difficult to verify the following 


lim Mx{{x,y)-, K) = Moo{{x,y)-, K) = < 
A —>-+00 


X' 


2 

2 

0 , 


— , if 0 < X < c, x = y, 


if x>c, X = y , 
if x^y, 


(5.14) 


with uniform convergence if y = x and x > 0. 

Despite its simplicity, this basic model elucidates the behaviour of Mx{{x,y); Kf) for a set Kg 
with a stair like boundary profile as, for instance, the one displayed in Figure [^a) . 

For such a set Kg, it is not difficult to verify that for (x, y) G let i G Z such that |(x — ic) — 
(y + ic)\ < c, then 

Mx{{x, y); K^) = Mx{{x -ic,y + ic); K'^) (5.15) 

with Mx{{x,y); K^) corresponding to the one step boundary domain perturbation discussed at the 
beginning of this example. Figure g6) contains the graph of Mx{{x,y); Kg), whereas Figure^c) 
shows its support, displayed together with the set Kg. The height of the ridges along Mk^ depends 
only on the gap size c, in particular, it is proportional to e?. It follows, therefore, that by setting the 
threshold larger than (?, the corresponding suplevel set of Mx{{x,y); Kg) will filter out all minor 
branches of Mk^, generated by the step-stair like boundary. 
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(a) 



1-5 


1.0 


0.5 


-0-5 


-1.0 l-l- ............................ i- 

-1-0 -0.5 0.0 0.5 1-0 1.5 2,0 

(b) 


Figure 8: Multiscale medial axis map of the set K displayed in Figure prototype of boundary 
perturbation as occurring in digitized images, (a) Graph of Mx{{x,y); K^) for A = 10; {h) Support 
of Mx{{x,y); K^) displayed together with the set K. 



Figure 9: (a) Geometric model of a prototype of staircase type perturbation of the boundary 

domain Kg. (b) Graph of Mx{{x,y); Kg) for A = 10; (c) Support of Mx{{x,y)‘, Kg) displayed 
together with the set K. 
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(a) (6) 


Figure 10: Multiscale medial axis map as obtained by the numerical implementation of 

Mx{{x,y)‘, K): (a) Support of the multiscale medial axis map Mx{{x,y)‘, K) for the digital im¬ 
age of a maple leaf, for A = 10. All the fine branches generated by the steps on the boundary are 
displayed.; (6) Suplevel set of Ma((x, y); K) corresponding to a threshold equal to one displaying 
only stable parts of the medial axis. 


As an application of these concepts, we show in Figure\T^ the results of the numerical realization 
of Mx{{x,y)] K) for the digital image of a maple leaf, where we can note the effects just discussed. 
In particular, Figure\Wla) depicts the support of Mx{{x,y)\ K) with the display of all fine branches 
created by the step-like irregularities of the boundary domain, whereas Figure \W{b) shows the su¬ 
plevel set of Mx{{x,y); K) corresponding to a threshold equal to one which singles out only the 
neighbourhood of stable parts of Mk- 


6 Proofs of Main Results 


Proof of Lemma 2.5: The existence of an affine function i such that ( |2.4[ )(i) and (2.4)(m) hold 
is well known (see e.g. m Remark 2.1] or [36]). The claim that 2 < A: < n -|- 1 comes from 
the Caratheodory Theorem [ITl Cor. 17.1.5 ] and the fact that co[/](0) < /(O). Also it is easy 
to see that xfs can be made distinct and Xi 0. For the proof of (2.4)(Hi), observe that / 
is upper semidifferentiable, i < f and i(xi) = f(xi). By [TS] Cor. 2.5], it thus follows that / 
is differentiable at xi and Df(xi) = Di(xi) = 2a. The proof of (2.4)(iv) is obtained from the 
definition of the convex envelope. For ( |2.4[ )(u), we have, by definition of the convex envelope, 
that I < co[/] and l'(O) = co[/](0). Since by hypothesis, co[/] G we can conclude that 

2a = Dl({)) = T»co[/](0), again using [HI Cor. 2.5]. □ 


Proof of Proposition 


2.13 


It is known (SH Theorem 3.1] that C')^(dist^(-; K)) G so 


we only need to improve the estimate of the Lipschitz constant obtained in |59[ pag. 755], namely 
8 -|- lOA. From the definition of the lower transform, we have 


C'i(dist^(-; K)){x + y)- C'i(dist^(-; K))(x) - (T>C'i(dist2(-; K))(x)) ■ y > -Alyl^ (6.1) 

for all x, y € M”. Now dist^(-; K) is a 2-semiconcave function |19[ Prop. 2.2.2], that is, x i—)• 
— dist^(a:, K) := g(x) is a convex function. So if we let fx(x) := (A -|- l)lxj^ — g(x), then by 
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Lemma 2.12, we have 


co[fx]{x + y) - co[fx]{x) - {D co[fx]{x)) ■ y < (A + l)|y|^ 
for all X, ?/ G M"'. We also have, for x G M”, that 

C'j^(dist^(-; K)){x) = co[dist^(-; iL) + A| • |^](x) — A|x|^ = co[/A(-)](ic) — A|x|^ . 


( 6 . 2 ) 


By (6.2), we obtain 

Ci(dist2(.; K))ix + y)- K)){x) - {DCi{dist\-, K))ix)) • y 

= co[fx]{x + y) -co[/a](x) - (L>co[/a](x)) ■ y - X\y\‘^ < (A + l)|yp - X\y\^ = |y|^, 


(6.3) 


Combining (6.1) and (6.3), we have 

-A|yp < C'i(dist2(-; K)){x + y) - C'i(dist2(-; K)){x) - {DC{{dlst‘^{-, K))(x)) ■ y < \y\^ . 

Thus X t-)- C'^(dist^(-; K)){x) is both 2-semiconcave and 2A-semiconvex. By [121 Corollary 3.3.8], 
we therefore conclude that C'^(dist^(-; K)) G C'^’^(M”'), and the Lipschitz constant of the gradient 
DC[{dist‘^{-; K)) is not greater than 2max{l, A}. □ 


Proof of Proposition 2.14: Let A > 0. Without loss of generality, we may assume that xq = 0. 
We consider two different cases, depending on the values of C'^(dist^(-; K)){0) and dist^(0; K). 

Case (i): C'(^(dist^(-; iil))(0) < dist^(0; K). 

In this case, by definition of the lower transform and the convex envelope, we have 

Ci(dist2(.; K)){0) = co[dist2(.; K) + X\ ■ |2](0). 

Since the function dist^(-; K) + A| • p is also continuous , up per semidifferentiable, coercive and 
co[dist^(-; K) X\ ■ p] G , the assumptions of Lemma 2.5 are satisfied. Let the affine function 
i and xi,..., Xfc G M”, Ai > 0,..., Afc > 0 with 2<A:<n-|-lbeas given by Lemma 2^, satisfying 
(2.4)(z) to (2.4)(x). We have, by (2.4 )(m) and (2.4)(iM), that 

2a • Xj + 6 = dist^(xi; iL)-|-A|xip, 2a = iA(dist^(-; iL) + A| • p)(xi). (6.4) 

Therefore x i—)• dist^(x; K) is differentiable at x* for i = 1, 2,...,/c. By [371 Lemma 8.5.12], we 
see that K{xi) := {yi} consists of a single element yi G K, so that dist^(xi; K) = \xi — ?/ip and 
Zldist^(-; K){xi) = 2(xj — yi). Thus ( |6.4[ ) reduces in this case to 

2a • Xj + 6 = |xi - yip + A|xip, 2a = 2(xi - yi) + 2Xxi. 


'^l\ 5 


By (6.5)2 and knowing that Yli=i ^ and X)i=i — 0) 

k k 


O = ^ Aif (Xi - yi) + Xxij = - ^ Ai^i, 


(6.5) 


( 6 . 6 ) 


i=l 


2=1 
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whereas by (6.5)i, (6.6) and the strict convexity of | • f, we obtain 
k k 

b = ^ Ai(2a • Xi + 6) = ^ Xi{\xi - yi\^ + A|xip) 


i=l 


2=1 

k 


> 


'^Xi{xi - Vi) 


2 = 1 


i=l 


= a 


Hence \a\^ < b, that is, 


(6.7) 


K)){0) 


<Ci(dist2(.; K)m. 


Case (ii): C^(dist^(-; K))(0) = dist^(0; K). 

In this case, we have 

co[dist^(-; K) + X\ ■ P](0) = C'^(dist^(-; i^))(0) = dist^(0; K) + A|0p, 


and 

co[dist^(-; K) + A| • |^](x) < dist^(x; K) + A|xp , 

for X G M"". Since co[dist^(-; K) + X\ ■ p] G is convex and dist^(-; K) + A| • p is upper- 
semidifferentiable, it follows from m Corollary 2.5] that D(dist^(x; K) + Ajxp) exists at 0, and 

Z)(dist2(-; K) + X\- |2)(0) = Z)co[dist2(.; i^) + A| • p](0). 

Again by m Lemma 8.5.12], K{0) = {uq} with uq G K the unique point that realizes the distance 
of 0 to K. So 

lZ)co[dist2(.; K) + X\- 12](0)12 = \D{dist\-, K) + X\ ■ l2)(0)l2 

= 1 — 2yol^ = 4dist^(0; K) = 4co[dist^(-; iL) + A] • 1^](0). 

□ 


Proof of Theorem 3.3^ Let A > 0. Without loss of generality, we may assume that xq = 0. We 
prove our result by establishing the contrapositive, and therefore suppose that 


Ci(dist2(-; K))i0) < dist2(0; K) 


(6.9) 


and seek to prove that 

Adist(0; Mk) < dist(0; K). 

As in the proof of Proposition 2.14^ all the assumptions of Lemma 2.5 are met for the function 
/(•) = dist^(-;iL) + A] • j^. Hence there exist an affine function £, points xi,...,Xk G M”, and 
Ai > 0,..., Afc > 0, 2 < A: < n + 1, that satisfy (2.4)(i) to (2.4)(n), which ensures that (6.5) holds 
with {yi} = K{xi)}, i G {1,..., k}. From (6.5 ) 2 , we have 


Xi = 


Vi + a 

1 + A 


( 6 . 10 ) 


32 






















which when substituted into (6.5)i yields 


2a ■ 


Vi + a 
1 + A 


+ b = 


Vi + a 


1 + A 


Vi 


+ A 


Vi + a 


1 + A 


A simple manipulation of this equation in m then gives 


Vi 


A 


2 (1 + A)|ap (1 + A)6 


A2 




and (6.10) and (6.11) together imply 


Xi - 


(1 + A)^ 


i = 1,2,, 


( 6 . 11 ) 


( 6 . 12 ) 


Now from ( |6.11 ), it follows that the points yi,y 2 , ■ ■ ■ ,yk G AT lie on the sphere S{a/X; c) = {?/ G 
M"", \y — a/A| = c}, and since Xi ^ Xj for i ^ j, we also have that yi ^ yj for i ^ j, by (6.10). We 
show next that the open ball B(a/X; c) does not intersect K, and hence a/A G Mk, the medial 
axis of K. We prove this claim by contradiction. Suppose y* G K (1 B{a/X] c), and define 


Then we have, from (2.4)(i), that 


X = 


2a-x* + b = £(x*} < dist2(x*, K) + A|x*|^ < |x* -y*\^ + A|x 


y* + a 
1 + A 


* 12 


.* 12 


12 


(6.13) 


(6.14) 


and by replacing (6.13) into (6.14), it follows that 

. a ' 


y 


X 


> c^ 


which contradicts the assumption that y* G B{a/X] c). Hence 


H(-;c) nK = 0, 


(6.15) 


and thus a/A G Mk- By the strict convexity of | • p and the fact that AiX, = 0, we then have, 

(1 + A)|ap (1 + X)b 


from (6.12), (6.7) and Proposition 2.14, that 

k 


disP^(0; Mk) < 


< 


i=l 


Xi - 


2 1 


(A+ 1)2 (A+ 1)2 

lap b 

+ 


A2 


' + 


A 


and hence 

This proves that if 0 ^ Mk and 


A2(l + A) A(l + A) A2 

Adist(0; Mk) < dist(0; K). 

^ ^ dist(0; K) 


b dist^(0; K) 

< T7T < 


(6.16) 


A2 


dist(0; Mk) ’ 

then C'(^(dist^(-; K)){0) > dist^(0; K). Since we always have that C'^(dist^(-; Ar))(0) < dist^(0; K), 
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it can be concluded that 

Ci(dist2(-; K))(0) = dist2(0; K), 

which completes the proof. □ 


Remark 6.1. Recall that in Remark 


3.1(b), we noted that if xq is a critical point o/C'3^(dist"^(-; K)), 


then Xq G co[iir(xo)]. Translating xq to 0, we can now see that i/C^(dist^(-; iir))(0) < dist^(0; iiT), 


this follows from ( |6.6| ), (6.11) and (6.15), since Lemma 2.5(v) implies thata = 0 if K )){0) 

0, whereas z/C^(dist^(-; K)){0) = dist^(0; ii'), the arguments in the proof of Proposition 2.14(ii) 
yield that 0 ^ K z/L>C'j^(dist^(-; iir))(0) = 0, thus clearly 0 G co[iir(0)] in this case also. 


Proof of Corollary 3.5| Note first that (3.7) and 


together yield that 


f| spprt(MA(-;K)) C f| Vx,k = Mk. 

0<A<+oo 0<A<+oo 

On the other hand, suppose xq G is such that M\{xq-, K) = 0. Then we have that dist^(xo; K) = 
C'j^(dist^(-; ii'))(xo), so as argued in the proof of Proposition 2.14 (ii), it follows that xq 0 Mk. Thus 
or all A > 0, Mk C {x G : Mx{x] K) >0}, which implies that 

Mk C. Pi spprtMA(-;7f). 

0<A<+oo 


□ 


Proof of Corollary |3.6[ We only need to consider the case where C'j^(dist^(-; K)){xo) < dist^(xo; K), 
since otherwise the claim is clearly true. Without loss of generality, assume that xq = 0. As in 
the proof of Theorem |3.3[ the assumptions of Lemma 2.5 are satisfied for the function /(•) = 
dist^(-; iL) + A| • p, so there exist an affine function £(x) = 2a • x + 6, points xi,..., Xfc G M”, and 
Ai > 0,..., Afc > 0, 2 < A: < n + 1, that satisfy (2.4)(z) to (2.4)(z;), so that (6.12) holds. Moreover, 
from Proposition 2.14 and the assumption that C'pdist^(-; K)){0) < dist^(0; K), it follows that 


< b < dist^(0; K). Hence, for each for z = 1, 2,..., A:, we have 


XU - 


|a| 

T 


< 


a 

X.-- 


1 + A 


< 


Vb 


and so 


|a| x/b 


2x/b ^ 2dist(0; K) 


A A 

Thus Xi G R(0; r-(O)) for z = 1, 2,..., A:, where r(0) = 2dist(0; K)/\. The conclusion then follows 
from the definitions of the convex envelope and of C'pdist^(-; iL))(0). 


□ 


Proof of Corollary 3.8: Since Mk C Vx^k, h follows from Theorem 3.3 that dist^(x;iL) = 
Cj^(dist^(-; K)){x) for all x in the set \ Vx^k, which is an open set because Vx^k is closed. Thus 
the result is immediate from Proposition 2.13 □ 


Proof of Proposition 3.11 Since dLl C then 

dist^(y,H'^) < dist^(z/,5H) for y G 


(6.17) 


34 































and by the ordering property of the lower transform, 


C'^(dist^(-; Q‘^))(x) < C^(dist^(-; dfi))(x) for x £ Q. (6.18) 

We want now to prove that the equality actually holds for x £ Q. We will show this by a contra¬ 
diction argument. Though the equality dist^(y; 11^) = dist^(y; 50) holds for all y £ fl, we cannot 
straightforwardly deduce the equality of the lower transforms in 0. Assume therefore, that at some 
point X G O we have that 


C(^(dist^(-; 0'^))(x) < C'^(dist^(-; 50))(x) 


(6.19) 


By the translation invariance of the distance and of the lower transform [641 Proposition 2.10], we 
can assume, without loss of generality, that x = 0 G O, so that (|6.19 ) becomes 


Ci(dist2(-; O‘=))(0) < Ci(dist2(-; 5O))(0), 
which is then equivalent to state that 

co[dist^(-; O'^) + A| • P](0) < co[dist^(-; 50) -|- A| • p](0). 


( 6 . 20 ) 


( 6 . 21 ) 


Since the function dist^(x,50) + A|xp is coercive and is continuous, by Proposition 2.3 Ai) there 
exists an affine function i(x) such that 


and 


i{x) < dist^(x; 50) + A|xp for all x G 


i{0) = co[dist^(-; 50) -|- A| • |^](0). 


Note that by Proposition 2.3 i), i{0) > 0. There must be a point yo £ (0)'^ such that 

i{yo) > dist^(yo; O^") + A|yoP = A|yoP- 


( 6 . 22 ) 

(6.23) 

(6.24) 


If we write our affine function as i{x) = a ■ x + b with b = £(0) given by (6.23) and a £ M”, then 

a-yo + b> A|yoP 


(6.24|) reads as 


(6.25) 


Since for x G 50, 


£{x) < dist^(x; 50) + A|xp = dist^(x; O'^) + A|xp , (6.26) 

and for t > 0 small enough, tyo £ O, it follows that there exists 0 < to < 1 for which toyo £ O'^ and 

£{toyo) = dist^(toyo; O'") X\toyo\'^ = A|toyo|^ • (6.27) 

This implies that 


too • yo + ^ = A|toyo| , 


that is. 


a-yo = toA|yo|^ - Vfo • 

If we substitute (6.29) into ( 6.25| ), we obtain 

toA|yo|^ - b/to + b> A|yo|^ , 


(6.28) 

(6.29) 

(6.30) 
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that is 


(6.31) 


b < -toX\yof < 0, 

which contradicts the fact that b > 0. Thus 

i{x) < dist^(x; 0'^) + A|x|^ for all x G M"', 

hence, 

co[dist^(-; + A| • P](0) > co[dist^(-; 50) + A| • |^](0) 

thus 

Ci(dist2(-; O‘=))(0) > Ci(dist2(-; 5O))(0), 
which contradicts the initial assumption (6.19). 


(6.32) 

(6.33) 

(6.34) 

□ 


Proof of Corollary 3.13[ We only need to verify that nO C In fact, if x G 14 , 0 '= hlO, 

then Adist(x; Mqc) < dist(x;0'^) < diam(O). Thus x G 114,o so that 0 \ 1 Ta,o C 0 \ 14 , 0 '=- The 
conclusion then follows from Corollary |3.8[ □ 


Proof of Theorem 3.15^ Note first that clearly M\[x] K) > 0 for all x G M”". We now prove 
the positive lower bound for M\{x] K) when x G Mk- Let x G Mk and r(x) = dist(x; K). Since 
K{x) C K then 

dist^(y; K) < dist^(y; K{x)) for all ?/ G M"' , 
hence, by the ordering property of the lower transform. Proposition 2.9 


C'^(dist^(-; K)){y) < C'^(dist^(-; K{x))){y) for all y . 
Now, by Proposition |2. 18 we have for y = x that 


(6.35) 

(6.36) 


Ci(dist2(-; iP(x)))(x) = (1 + A)dist2(x; co[iL 2 ,A(a;)]) + ■ 


1 +A 


(6.37) 


By the Caratheodory’s theorem m, for every w = x + y & co[iL2,A(3;)], there are at most n + 1 
points X + yi/{l + A) G K 2 ^\{x) with x + ?/* G K{x), i.e. x + y* G dK and \yi\ = r(x), and 
Ai > 0,..., An+i > 0 such that ^ w = x + + ^)- Thus we have 


X — rc = 


n+l 


i=l 


Vi 


+ A 


(1 + A)^ 


n+l 


2=1 


(1 + A)^ 


n+l 


X — 


^Aj(x + y4 


2=1 


which yields 


dist^(x; co[iL 2 ,A(a:^)]) = , ^ „ dist^(x; co[iP(x)]). 

(i + A) 


(6.38) 


(6.39) 


By substituting (6.39) into (6.37) we have 


Ci(dist2(-; K{x))){x) = ^^dist2(x; co[iL(x)]) + ^^r^ix) 


(6.40) 


By comparing ( |6.36 ) and (6.40), we finally obtain 

Ma(x; iL) > (1 + A)(dist^(x; K) — C'^(dist^(x; K{x)))) = dist^(x; K) — dist^(x; co[iL(x)]) 
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To find now an upper bound to M\{x] K) that holds for all x G M”, note first that x € K, then 
Ma(x; K) = dist^(x; K) = 0. Suppose now that x ^ K. Then r(x) > 0 and B{x-, r(x)) C so 
K C B^{x]r{x)) and hence 

dist^(y; B^{x,r{x))) < dist^(y; K) for y G M” , (6-41) 

and thus, by the ordering property of the lower transform (Proposition |2.9[ ), 

C'i(dist^(-;-B‘'(x,r(x))))(y) < C'i(dist^(-; i^))(y) for y G M'" . (6.42) 


By Lemma 2.20, for |y| < r(x)/(l + A), after a simple translation of points and due to the invariance 
of the distance transform, we have 


which for y = x gives 


C'i(dist^(-; B^{x, r(x))))(y) = “ ^1^ 


Ci(dist2(-; 5"(x,r(x))))(x) = 


hence 


dist^(x; K) — C'^(dist^(-; B^{x,r{x)))){x) = dist^(x; K). 

1 + A 


(6.43) 

(6.44) 

(6.45) 


By comparing (6.42) and (6.45), we then conclude that 


dist^(x; K) — C'^(dist^(-; K)){x) < dist^(x; K) — C'j^(dist^(-; B^{x-, r(x))))(x) = dist^(x; K ), 

1 + A 


which completes the proof. 


□ 


Proof of Proposition 3.18^ Let X G ilLjc, r(x) — dist(x 5 ^ 0, and denote by x^, X2 G 
B{x; r(x)) n K the points of K{x) that realize the separation angle Ox at the point x. Thus 
|x — xi| = |x — X 2 I = r(x). Since {xi,X 2 } C K(x), then xi + t{x 2 — xi) G co[iL(x)] for t G [0,1], 
which for t = 1/2 gives 


+ X2. f Ox \ . 

-z-) = cos — r(x). 


Thus 


and hence 


dist(x; co[iL(x)]) < dist(x; 


dist^(x; co[iL(x)]) < dist^(x; K) cos^ ( ^ 


dist^(x; K) — dist^(x; co[iL(x)]) > dist^(x; K) (1 — cos^ ( % )) = dist^(x; K) sin^ ( ^ ) > (6.48) 


2 / 


(6.46) 

(6.47) 


as required. 


□ 


Proof of Proposition 3.20: If xq G K, clearly Mx{xq', K) = 0 for all A > 0. So we may assume 
that xo ^ K. Since xq ^ Mk, dist^(-; K) is differentiable at xq [371 Lemma 8.5.12]. Therefore for 
every e > 0, there exists (5 > 0 such that 


|dist^(xo + y; iL) - dist^(xo; K) - Ddist^(xo; K) ■ y\ < ejyj 
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3.6 


we have 


for y G i?(0; S). Now by the locality property Corollary 

Ci(dist2(-; K))(xo) = K) + A|(-) - xoP](xo) 

where r(xo) = 2dist(xo; K)/X > 0. Thus for A > 0 sufficien tly large, r(xo) < (1. Since x i—>■ 
dist^(x; K) + A|x — xqP is continuous and coercive, by Lemma 2.5 and Corollary 3.6, there exist 
xi,...,Xk G B{xo; r{xo)) and Ai > 0,..., A^ > 0 such that Yli=i = 1) Y^i=i and 


co.B(xo;r(xo))[dist^(-; + Mi') “ a^o|^](a:o) = Yj^^ {dist^{xi; K) + A|xi - xo|^) 

i=l 

k 

> ^ Aj (dist^(a;o; K) + iAdist^(xo; K) ■ (x* - xq) - e|xi - xo| + A|xj - xop) 


i=l 


> dist^(xo; K) — 


4A 


as |xi — xo| < r(xo) < 5. Here we have also used the facts that AiiAdist^(xo; K) ■ (xj — xq) = 0 
and that Xt^ — et > —e^/(4A) for t G M. Since we also have C'^(dist^(-; K)){xo) < dist^(xo; K), we 
have 

0<Mx{xo-,K)< + 


Thus 


4A 


0 < liminfMA(xo; K) < hmsupMA(xo; K) < — . 

A— >+CXD A— >-+00 4: 


Since e > 0 is arbitrary, the conclusion follows. 


□ 


Proof of Theorem 3.23: We only consider the case xq G Mk- Again without loss of generality, 
we may assume that xq = 0 G Mk- Let Kq = Ar(0) and dist(0; K) = dist(0; Kq) = tq > 0. Since 
Ko C K, we have dist^(x; K) < dist^(x; Kq) for x G M"', so that 


C'i(dist^(-; K)){x) < C'i(dist^(-; Ko)){x) = 


^^0 

1 + A 


+ (1 + A)dist^(x; co[A'o/(l + -^)]) ~ (6.49) 


for X G M"'. Therefore 

Ma(0; K) > Ma(0; Kq) = (1 + A)(dist2(0; Kq) - Ko)){0)) 

= rg - dist2(0; co[Ko]) = M(0; Ko). 


(6.50) 


3.6 


Next we establish lower bounds for C'(^(dist^(-; K)){x) using the locality property from Corollary 
For 0 < e < ro sufficiently small, let Kq-^ C ^(O; xq) be the closed e-neighbourhood of 


Kq on the sphere, defined using the geodesic distance p(x, y) on 5(0; xq), that is iLo;e = {u ^ 
5'(0; xo), p(y; Kq) < e}, where p(y; Kq) = inf{p(y,x), x G Kq}. 

The aim of the following technical construction is to show that for x in a small neighbourhood 
of 0, dist^(x; Kq, 2 e) is a lower bound for dist^(x; K). For <5 > 0, define the closed neighbourhood 
Kq^ = {(xq + t)y/\y\, y G Kq^^, 0 < t < <5} and note that ATg ^ is clearly a compact set. Then it can 
easily be proved, using a contradiction argument, that for every 0 < e < 1, there exists 0 < J 
such that K n H(0; xq + (5) C K^^. Define also another compact set by 


= Kl, U iLo,2., 
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where i^o, 2 e will be used to ‘shadow’ and the unbounded closed set 

= ro + h). 

Clearly, K C so that dist^(x; K) > dist^(x; W^^-s) for all x G M"’. 
We claim that there exists rj > 0 snfficiently small such that 


dist^(x; W^-s) = dist^(x; Ko, 2 e) 


(6.51) 


for X G B{0-, rf). We postpone the proof of (6.51) to the end and proceed first to assume that (6.51) 
holds. Then for A > 0 sufficiently large, we have 


2dist(0; W^^s) 2dist(0; il"o, 2 e) 2ro 


< r]. 


A A 

By the locality property (Corollary |3.6[ ), we have 

Ci(dist2(.; W,,5))(0) - (dist2(.; + A| • |2)(0) 

= co:B(0;,,)[dist2(.; W,^s) + A| • p)](0) 

= COB(0;r,)[dX(-; Ko.2e) + A| • |2](0) 

= co[dist^(-; KQ, 2 e) + A| • |^](0) 

= Ci(dist2(.;ii:o,2.))(0) 

X,2e 


\r^ 

-—^ + (1 + A)dist^ ( 0; co 
1 + A 


1 + A 


where we have used (6.51) and Proposition 2.19 Thus we obtain 


Ci(dist2(.; K)m > Ci(dist2(.; W,,,))(0) = 


1 + A 


(6.52) 


+ (1 + A)dist^(0; co[iA’o, 2 e/(l + A)]). 


As dist^(0; K) = dist^(0; = dist^(0; A'o, 2 e) = fQ, we then have 

Ma(0; K) <rl- dist^(0; co[A:o, 2 e]) ■ 
Therefore for sufficiently large A > 0, 


Moo(0; Kq) < M\{0; K) < Moo(0; Ko^2e)- (6.53) 

Passing to the limit A —)• +oo then gives that for each fixed e > 0 small, 

Moo(0; Ko) < liminf Ma(0; K) < limsupMA(0; K) < Moo(0; Aro, 2 e) ■ 

A^+oo A^+cxD 

Since Kq is compact and Kq ,2e Kq as e —)• 0 under the Hausdorff distance in we also 

have that co[A'o, 2 e] —^ co[A'o] as e —>■ 0 under the Hausdorff distance in M"". Thus as H i—>■ 
dist^(0; V) is continuous under the Hausdorff distance for compact sets V C [6], it follows 
that hme_^o+ dist^(0; co[iiro, 2 e]) = dist^(0; co[iCo]), and hence lime_^o+ Afoo(0; KQ^ 2 e) = Afoo(0; Kq). 
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Hence limA-s.+oo-^ a(0; K) exists, and 


lim Ma(0; K) = M^{0; Ko) = M^{0- K). 

A—>-+oo 


It remains to prove (6.51). First note that when 0 < rj < 6/2, 


dist(x; 5(0; tq + (J)) > dist(x; iiro, 2 e) 


(6.54) 


because dist(x; 5(0; ro+d)) = ro+(5—|x| > ro+(5 —t/ and dist(x; i^o, 2 e) < dist(0; ii'o, 2 e)+|a^| < ro+r], 


so that (6.54) holds if + t] < vq + 6 — rj, which is equivalent to 2r] < 6. 

Now we show that dist(x; = dist(x; ^Co, 2 e)- Given any point zq = (to + 'f'o)yo/\vo\ £ 
Vf:^s \ KQ^ 2 e, with 0 < to < 6 and yo G Kq^^, we observe that a necessary condition for some 
X G H(0; rj) to reach the distance to 14^5 at zo, that is, dist(x; V^^s) = \x — zo\, is that the line passing 
through Zo and x does not intersect iFo, 2 e- Notice that for the point yo G Ko^e, the e-neighbourhood 
of yo in 5(0; tq) under the geodesic distance p, given by := {tc G 5(0; ro), p{yo,w) < e} is 
contained in iFo, 2 e- Therefore if we draw a line passing through zo and the relative boundary of 
5yp_e in 5(0; tq) and we can show that the distance between the line and the origin 0 is bounded 
below by a positive constant uniformly with respect to yo G Ko,e and 0 < to < 6, then we can find 
0 < r] < 5/2, such that dist(a:; 14 ^ 5 ) = dist(x; Kope) for x G B{0; p). 

Due to the symmetry of Euclidean balls and spheres, we only need to consider the case in with 
2/0 = (^ 0 ) 0), z = (r-Q+t, 0), where 0 < t < 6 and = {(^0 cos 6, ro sin0), —e/rg < 0 < e/ro}- The 
distance between the line L passing through z and the boundary point (ro cos(e/ro), rosin(e/ro)) 
and the origin (0,0) is attained at a point of the form (s,tt(s)), where 


u{s) = 


ro sin(e/ro)(ro + t- s) 


ro + t — r cos(e/ro) ’ 
so that the squared-distance between (0,0) and a point (s,tt(s)) in L is 

rg sin2(e/ro)(ro -ht - s)^ 


+ 


(ro -I- t - r cos(e/ro))2 


with the minimum point at 

_ _ rg sin^(e/ro)(ro +t) _ 

r^ sin^(e/ro) -|- (ro -|- t - r cos(e/ro ))2 

The distance between (0, 0) and L is 



+ u'^iso) > |so| > 


rQ sin^(e/ro) 
2 (ro + <5)2 


:= ?7o > 0 . 


Therefore if we choose 0 < p < min{(5/2, po}, we have, for all x G B{0; p), that dist(x; 14 ^, 5 ) = 
dist(a:; Kope), and hence 

dist(x; W^^s) = dist(x; Kop^) 

for all X G B{0; p). □ 


Proof of Proposition 4.1: This follows from the dehnition of M\{{x,y)', K) and Lemma 2.21 


□ 
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Proof of Theorem 4.3: Let fj. = dist'^(iL; L) with /i finite since K and L are compact sets. By 
Definition (2.23) for Hausdorff distance, we have for x G M” 

|dist(x; K) — dist(x; L)\ < fi, (6.55) 

hence, 

dist^(x; K) < + 2/rdist(x; L) + dist^(x; L) 

< /i^ + //(I + dist^(x; L)) + dist^(x; L) = /i(l + /r) + (1 + /i)dist^(x; L). 


After adding A|xp to both sides and taking the convex envelope, we find 


co[dist^(-; AT) + A| • |^](x) < ^(1 + /r) + (1 + /r) co[dist^(-; L) + 


1 + /i 


^](x) 


which yields 


Ci(dist2(-; K)){x) < ^(1 + //) + (! + ^)C^/(i+^)(dist2(-; L)){x) 


(6.56) 


Since 


<^A/(i+M)(dist^(-; L)){x) < Ci(dist2(-; L))(rE) and Cj^/(i+^)(dist2(-; L)){x) < dist2(x, L) 


we obtain, from (6.56), after using (6.55), that 

C'^(dist^(-; K)){x) < /i(l + ^) + C^(dist^(-; L)){x) + ^dist^(a:; L) 

< C'^(dist^(-; L)){x) + /i(l + /r) + fi{fL + dist(x; K) f). 
With a similar argument, we find that 

C'^(dist^(-; L)){x) < ^(1 + /r) + C'^(dist^(-; K)){x) + ^dist^(a:; K) 

< C'^(dist^(-; K)){x) + ^(1 + //) + + dist(x; K))"^). 


(6.57) 


(6.58) 


By comparing (6.57) and (6.58) we therefore conclude that given a compact set K C M”, for any 
compact set L C W^, we have that, for any x G M”, 

C'^(dist^(-; K)){x) - C{{disi‘^{-; L)){x) <//(^(l + ^) + (/r + dist(x; AT))^)^ , (6.59) 


which proves (4.2). To show ( |4.3| ), observe that after using (6.55) we have for any x G M” 
|dist^(x; AT) — dist^(x; L)\ < |dist(x; K) — dist(x; A)||dist(x; AT) + dist(x; L)\ 

< /x(2dist(x; K) + /r), 


(6.60) 


and from the definition of the multiscale medial axis map and the triangle inequality, we obtain 
|Ma(x; K) — M\{x‘, L)\ < ^(1 + A) ^(dist(x; AT) + + 2dist(x; K) + 2^ + 1^ , (6.61) 


where we have taken into account (6.59) and (6.60). This concludes the proof. 


□ 
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since dist(x; dQ) < diam(r2) if x G 

□ 
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Proof of Corollary 


4.5 


This follows from Theorem 4.3 
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